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Though the dependence of near-field radiative transfer on the gap between two planar objects is 
well understood, that between curved objects is still unclear. We show, based on the analysis of the 
surface polariton mediated radiative transfer between two spheres of equal radii R and minimum 
gap d, that the near-field radiative transfer scales as R/d as d/R — » and as In (R/d) for larger 
values of d/R up to the far-field limit. We propose a modified form of the proximity approximation 
to predict near-field radiative transfer between curved objects from simulations of radiative transfer 
between planar surfaces. 
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For more than a century, Planck's theory of black- 
body radiation and the radiative transfer theory (RTT) 
have been successful at predicting radiative heat trans- 
fer between objects when all length scales involved are 
much larger than the characteristic thermal wavelength 
At given by Wien's displacement law (At ~ 1.27fk3//c#T, 
where 2irh is the Planck's constant, c is the speed of 
light in vacuum, ks is the Boltzmann constant and T is 
the absolute temperature). However, at shorter length 
scales, radiative transfer between two objects can exceed 
the predictions of Planck's theory due to electromagnetic 
near-field effects. Recent precision experiments^"^ for 
measuring radiative heat transfer between a silica mi- 
crosphere and a planar silica substrate have begun to 
shed new light on the enhancement of radiative trans- 
fer at nanoscale gaps due to surface phonon-polaritons. 
While near-field enhancement between parallel surfaces 
is well understood theoretically,'^' that between curved 
surfaces is still unclear. Though near-field radiative ex- 
change between a nanosphere and a plane has been cal- 
culated under different approximations p03 rigorous nu- 
merical computation of near-field radiative transfer be- 
tween a sphere of arbitrary diameter and a planar surface 
has not yet been possible due to computational difficul- 
ties. In order to understand the effect of curvature on 
enhancement of radiative transfer at nanoscale gaps, we 
investigate the near-field radiative transfer between two 
spheres of equal radii by rigorous simulations based on 
the dyadic Green's function technique^ and fluctuation- 
dissipation theoremfm 

The problem of computing fluctuation-induced inter- 
action between objects with curved surfaces is also en- 
countered in literature on van der Waals and Casimir 
forces. A frequently used strategy to obtain forces be- 
tween spherical surfaces, when the force between two 
parallel surfaces is known as a function of separation, 
is to use the so-called proximity approximation (or Der- 
jaguin approximation). 12 While a rigorous proof for the 
proximity approximation for forces is still lacking,^ it's 
accuracy in predicting Casimir force has been well in- 
vestigated theoreticalf^^tti^] and is expected to be most 
accurate when R ^> d, where R is a characteristic ra- 



dius of the objects involved, and d is the minimum gap 
between them. Experimental measurements of Casimir 
forc d 17 * 1 ^ have confirmed the applicability of the proxim- 
ity approximation when R^> d. 

A similar technique^ has been used to compute radia- 
tive transfer. Measurements from the Chen grou 
tween a silica microsphere and a silica substrate in the 
range 30 nm to 10 fim seemed to show better agreement 
with theoretical predictions of heat transfer between two 
spheres^ and did not agree with the proximity approx- 
imation. Rousseau et al.pl based on their measurements 
between a silica microsphere and a silica substrate, con- 
cluded that near-field radiative transfer agreed with the 
proximity approximation in the range 30 nm to 2.5 /am. 
There are no experiments between two silica spheres re- 
ported in literature. Though the phenomena of van der 
Waals force (including Casimir force) and near-field ra- 
diative transfer are fluctuation-induced, there are im- 
portant differences. Radiative transfer has contributions 
from the infra-red (IR) portion of the electromagnetic 
spectrum whereas forces have larger contributions from 
the visible and higher frequencies. Dispersion forces obey 
a power law behavior and decay rapidly to zero as gap 
between the interacting bodies increases, while radiative 
transfer has a finite value due to propagating waves at 
large gaps too. Because of these differences, it is not 
clear whether the proximity approximation, as it is used 
to compute dispersion forces, can be used to predict near- 
field radiative heat transfer between spherical surfaces. 

The configuration of the two spheres between which 
radiative transfer is to be calculated is shown in [T] 
(top right corner). Radiative heat transfer between the 
spheres is calculated using Rytov's theory of fluctuational 
electrodynamics.^ The Fourier component of the fluctu- 
ating electric field E{r\,uj) and magnetic field H(r\,ui) 
at any point r\ is given by^ 

E(n,Lj) = iujfio d 3 rG e (ri,r,uj) • J(r,u), (1) 
Jv 

H(n,u) = [ d 3 r^ h (r u r^)-J(r^), (2) 
Jv 

where G e (r , i,r,o;) and G^(ri, r, ui) are the dyadic 
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Green's functions for the electric and magnetic fields due 

to a point source at r and are related by G^(ri,r,o;) = 

V x G e (ri, r, u;), J{r,uj) is the Fourier component of 
the current density due to thermal fluctuations, and fi Q 
is the permeability of vacuum. The integration is per- 
formed over the entire volume V containing the source. 

The ensemble average Poynting vector S = \ Re{E x 
H*) depends on the cross spectral densities of the compo- 
nents of the fluctuating current source, which are related 
to temperature by the fluctuation-dissipation theorem. 
Further details on the analysis can be found in Ref. 19. 
Refinement of the numerical method has enabled us to 
probe lower gaps and suggest better models for near-field 
radiative transfer. To investigate the effects of surface 
phonon-polaritons on near-field radiative transfer, silica 
has been the material of choice in experiments for two 
reasons: (1) it can support surface phonon-polaritons in 
the frequency ranges from 0.055 to 0.07 eV and 0.114 to 
0.16 eV, and (2) silica microspheres are easily available. 
Hence the heat transfer has been computed, using silica 
as the material, for the frequency range 0.041 eV to 0.16 
eV. The optical properties of silica are taken from Ref. 
20. All numerical simulations have been conducted at 
300 K. 

The linearized thermal conductance G (WK _1 ) be- 
tween the two spheres is defined as: 



G= lim %™\ 

tht 2 T 1 -T 2 



(3) 



where Q (Xi, T2) is the rate of heat transfer between the 
two spheres at temperatures T\ and T 2 . Numerical val- 
ues of conductance are plotted as a function of gap to 
radius ratio d/R for different radii in [I] Only the points 
which deviate significantly from RTT have been included 
(see supplementary information, Section 1, for further 
details). In [I] for every radius, two regions can be ob- 
served: (1) a region where conductance varies logarith- 
mically (marked Region-A), and (2) a region where a 
deviation from logarithmic behavior is observed (marked 
Region-B). To gain a deeper insight into this behavior, 
we compare the spectral variation of the conductance at 
different gaps in the two regions. The gaps chosen are 
d/R = 0.01, 0.03 and 0.05 for R = 20 /mi (marked (a), 
(b) and (c) in [TJ. The gaps corresponding to d/R = 
0.03, 0.05 fall in Region-A of [TJ while d/R = 0.01 falls 
in Region-B. From the spectral variation plotted in |2j 
we note that most of the increase in the heat transfer 
for d/R = 0.01 is due to the contributions from surface 
phonon-polaritons alone. Hence we conclude that the 
contribution from surface phonon-polaritons to the con- 
ductance starts to become significant in Region-B. 

We also analyzed the contributions of conductance 
from the resonant and nonresonant frequencies sepa- 
rately. The results are shown in [3] Remarkably, the 
analysis for the spectral conductance at a resonant fre- 
quency (0.061 eV) suggests that at gaps d/R < 0.01, the 
conductance is dependent only on the ratio d/R and is 



CD 
O 

B 

o 

=3 
"D 

o 
O 



1 11 1 

T 


11 1 

R^. 






Region - b\ 




• \ 

\ * 


V(a) Nz Re 9' on " A 


t\ 




■ R = 5 nm V 




♦ /?=10um\ 




• R = 15 urn 




A # = 20 jam 




▼ R = 25 urn 





10~ 3 10 " 1 ' 
Nondimensional gap (d/R) 



v2 



FIG. 1: Conductance between the two spheres shown in the 
top right corner as a function of d/R for different radii. The 
open circles denote the conductance values which show a loga- 
rithmic variation with gap (marked Region-A) and the closed 
circles denote the conductance values which show a deviation 
from logarithmic behavior (marked Region-B). The spectral 
variation of the conductance at gaps marked (a), (b) and (c) 
for R = 20 /im spheres are shown in [2] 
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FIG. 2: The spectral variation of conductance for R = 20 fim 
for the different gaps (a), (b) and (c) marked in [l] The 
frequency regions marked "Resonant frequencies" ("Nonres- 
onant frequencies") are where surface phonon-polaritons are 
present (absent). 
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FIG. 3: For different radii, the spectral conductance at a 
resonant frequency (0.061 eV) as a function of d/R. The 
conductance values for all radii attain a slope of —1 at low 
gaps. Inset: The spectral conductance at a nonresonant fre- 
quency (0.1005 eV) as a function of d/R (axes labels remain 
the same). 



independent of the particular values of d and R. Further- 
more, the slope of the data points being w — 1 suggests 
a (R/d) behavior at such gaps. A similar analysis for a 
nonresonant frequency (0.1005 eV), shown in the inset of 
[3j suggests that for d/R < 0.01 the rate of change of spec- 
tral conductance with gap is significantly lower than that 
for resonant frequencies. Based on the behavior of reso- 
nant radiative transfer at small gaps and the observation 
of logarithmic behavior for larger gaps in[l] we have found 
that the numerical values of conductance can be modeled 
by a function of the form d (R/d) + C 2 In (R/d) + C 3 , 
where Ci, C2 and C3 are radius-dependent constants. 
Ci, C2 and Cs are obtained by minimizing the square 
of error between the function and the numerical values 
of conductance. The fitted curves are shown in [I] (For 
details on the variation of C±, C2 and C3 with R see 
supplementary material, Fig. 2). 

As has been mentioned earlier, conductance between 
curved surfaces can also be estimated from the known 
solutions for near-field radiative transfer between paral- 
lel surfaces using the proximity approximation. The heat 
transfer coefficient h(z) for two flat silica surfaces is plot- 
ted as a function of gap z inQa). h(z) can be split as 
follows: 



h(z) = h n f (z) + h Q 



(4) 



where is the contribution from propagating waves, 
which attains a constant value for z ^> At, and h n f con- 



tains contributions from all other effects. As can be seen 
from[4ja), h n f (z) — >• as z —> 00. 

Rousseau et alP applied the proximity approximation 
to compute the theoretical conductance between a sphere 
and a flat surface. Applying the same form of proximity 
approximation for finding the conductance between two 
spheres, we get: 



G(d, T) 



h(z)2i\r dr, 



(5) 



where z — d + 2R — 2\/R? — r 2 is the local gap be- 
tween the two spheres as shown in [4^b) and h(z) is the 
heat transfer coefficient between two parallel surfaces at 
that gap. From[4ja), we observe that h n f » for 
gaps z < 400 nm. When sizes of spheres are such that 
z(r) < 400 nm for < r < R, most of the contribu- 
tion to the conductance is from near-field effects. How- 
ever, for sizes of spheres currently used in experiments 
(R > 5 jam) there will be significant contribution from 
larger gaps where near-field effects are minimal. While 
[5] might provide a reasonable approximation to the con- 
tribution from near-field effects to the radiative trans- 
fer, it overestimates the contribution from propagating 
waves. In the large-gap limit, where the heat transfer 
coefficient attains a constant value, h^, [5] predicts con- 
ductance between the spheres to be i:R 2 h 00l irrespective 
of the gap and does not take into account the variation 
of view factoid between the two spheres with distance. 

In lieu of this we propose that the proximity approx- 
imation, in the form that is used to determine Casimir 
or van der Waals forces between spherical surfaces, be 
used to predict only the contribution to conductance from 
h n f. The contribution to the conductance from propa- 
gating waves is computed according to RTT by taking 
into account the variation of view factor with distance 
between the two spheres. This correction to the prox- 
imity approximation formulation is not necessary while 
calculating Casimir or van der Waals force, since they de- 
cay rapidly with distance (l/d A and 1/d 3 respectively). 
The modified form of proximity approximation to deter- 
mine the conductance is: 



G(d,T) = [ h nf (z)27rrdr + G c (d,T), 
Jo 



(6) 



where G c (d, T) can be approximated by the conductance 
value from RTT when diffraction effects are negligible. 
G c (d, T) for two objects of equal emissivity e and surface 
area A is given by^: 



G c (d,T) 



4(7 AT 3 



2(1 -e)/e+(l/ F 12 y 



(7) 



where F12 is the view factor between the two objects. 
Conductance values computed using |6] and [7| are m 
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FIG. 4: (a) Variation of h(z) and h n f(z) with gap between 
two flat silica surfaces, (b) Proximity approximation - the 
conductance between two spheres is calculated by summing 
the local contributions of the heat transfer coefficient between 
two parallel planes, (c) Comparison between the conductance 
values obtained numerically and using proximity approxima- 
tion. The open cirlces denote the numerical values while the 
closed circles denote the proximity approximation predictions 
using [6] Proximity approximation using [5] for R = 25 /xm has 
been included for comparison. 



greater agreement with the numerically computed val- 
ues than the prediction of [5] as shown in Qc). Since 



the conductance values have been computed for the fre- 
quency range 0.041 eV to 0.164 eV, G c (d,T) has been 
determined for this frequency range too (see supplemen- 
tary material, Section 1, for more details). Numerical 
values of the gap dependent view factor between the two 
spheres is taken from Ref. 22. 

While the modified form of the proximity approxima- 
tion can also be used to predict the conductance between 
a sphere and a flat surface, there are no rigorous numer- 
ical simulations of near-field radiative transfer between 
a sphere and a flat surface to compare them with. If 
this can be accomplished, the validity of this modified 
proximity approximation for predicting the conductance 
between sphere and a plane can be verified. Another 
configuration that can be used to test the predictive ca- 
pability of the modified proximity approximation is that 
between two parallel cylinders. Current experimental se- 
tups are not ideal for comparison either. Since the sphere 
is positioned close to the edge of the substrate due to ex- 
perimental constraints, deviation from RTT could occur 
due to surface phonon-polariton mediated near-field ef- 
fects as well as diffraction effects which need to be com- 
puted for each geometry separately. 

In summary, radiative heat transfer between two 
spheres has been analyzed in the near-field regime us- 
ing flucutational electrodynamics formalism. We have 
shown that it varies as R/d as d/R — ^ and as In (R/d) 
for larger values of d/R up to the far-field limit. We have 
also shown that the proximity theorem, in the form that 
is used to compute dispersive forces cannot be used to 
determine near-field heat transfer and a modification is 
needed to take into account the contributions from prop- 
agating waves. Further numerical simulations of near- 
field radiative transfer between curved surfaces are nec- 
essary to test the predictive capability of the modified 
proximity approximation in more general configurations 
than the one considered here. Such predictive capabilities 
will be useful in applications like heat-assisted magnetic 
recording^ to approximate the near-field radiative trans- 
fer between the curved surface of the near-field trans- 
ducer and the recording medium. 
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